Mon. Not. R. Astron. Soc. 000, [T]f5] (2012) Printed 19 July 2012 (MN lAT^K style file v2.2) 



On the Evolution of the Snow Line in Protoplanetary Discs 



Rebecca G. Martin and Mario Livio 

Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 



(N 

o 

(N 
00 



Oh; 

6 

o3 



> 

00 

o 

(N 



ABSTRACT 

We model the evolution of the snow line in a protoplanetary disc. If the magneto- 
rotational instability (MRI) drives turbulence throughout the disc, there is a unique 
snow line outside of which the disc is icy. The snow line moves closer to the star as the 
infall accretion rate drops. Because the snow line moves inside the radius of the Earth's 
orbit, the formation of our water-devoid planet is difficult with this model. However, 
protoplanetary discs are not likely to be sufficiently ionised to be fully turbulent. A 
dead zone at the mid-plane slows the flow of material through the disc and a steady 
state cannot be achieved. We therefore model the evolution of the snow line also 
in a time-dependent disc with a dead zone. As the mass is accumulating, the outer 
parts of the dead zone become self gravitating, heat the massive disc and thus the 
outer snow line does not come inside the radius of the Earth's orbit, contrary to 
the fully turbulent disc model. There is a second, inner icy region, within the dead 
zone, that moves inwards of the Earth's orbit after a time of about 10 6 yr. With this 
model there is sufficient time and mass in the disc for the Earth to form from water- 
devoid planetesimals at a radius of 1 AU. Furthermore, the additional inner icy region 
predicted by this model may allow for the formation of giant planets close to their 
host star without the need for much migration. 

Key words: accretion, accretion discs - protoplanetary discs - stars: pre-main- 
sequence - planets and satellites: formation - Earth 



1 INTRODUCTION 

Bodies in our solar system show a distribution of water 
abundance. The innermost terrestrial planets contain little 
water compared with the outer planets. Planet formation 
is thought to occur from planetesimals that formed in the 
solar nebula. Within the solar nebula, ice forms beyond a 
radius from the central star known as the snow line, iisnow- 
This is thought to play an important role in the composition 
of forming planets. The solid mass density outside the snow 
line is much higher because of water ice condensation. In 
fact, for a solar compos ition, water ice abu ndance is as high 
as silicate and iron (e.g. iPollack et ai1ll994h . Observations of 
the asteroid belt, located between Mars and Jupiter, suggest 
that the snow line is currently located within t his region. 
The outer asteroids are icy C-class objects (e.g. lAbe et all 
l200d : lMorbidelli et alj|2000h whereas the inner asteroid belt 
is largely devoid of water. This implies that when planetesi- 
mal formation occurred the snow line was located at around 
-Rsnow = 2.7 AU from the Sun. 

The snow line occurs at a te mperature, T Bnovi , 
that is in the range of 145 K i|Podolak fc Zuckerl 
l2004h to 170 K (|Havashil Il98l"l ), depending on the par- 
tial pressure of nebular water vapour. Its distance 
from the star is usually calculated in a fully turbu- 



lent steady state accretion disc (e.g. ISasselov fc Lecarl 
20001: iLecar et al.ll2006l; iKennedv. Kenvon fc Bromlevll2006l ; 



Kennedy fc Kenvonl 20081 ). The magneto-rotational insta- 
bility (MRI) is thought to drive the turbulence in the 
disc (jBalbus fc Hawlevl [l991) and thus transport angular 
moment um outwards allowing ac cretion on to the central 
star (e.g. lFromang fc Nel son 2006). Recent calculations have 
used a detailed disc structure and include a stellar radi- 
ation flux and viscous dissipation of the gas as the main 
heating sources. The disc is assumed to be in a quasi- 
stea dy state as the i nfall accretion rate decreases in tim e 
(e.g.llda fc Lm1l2005l : iGaraud fc Lm]|2007l : iMin et al.ll2ulll ). 
iDavisI l|2005l ) found that the radius of the snow line reaches 
a minimum of about 0.6 AU, insid e the current orbit o f 
the Earth and t his w as confirmed by IGaraud fc Linl l|2007h . 
IGaraud fc Lml (|2007l ) find that the snow line migrates in- 
wards as the accretion rate drops (down to an accretion rate 
of around M — 10~ 10 Mq yr -1 ) because the viscous dissipa- 
tion decreases. However, as the accretion rate drops below 
lO" lo M yr-\ the snow line migrates outwards again as 
the disc becomes optically thin and the temperature rises. 

The fully turbulent steady state model predicts a mini- 
mum snow line radius within the Earth's orbit. If planetesi- 
mals formed during this time, the Earth would have formed 
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from icy bodies. This appears to be contradictory with the 
current water cont ent on Earth that is very low at around 
0.023% by weig ht <|Lewisll2004h . For comparison, the outer 
solar system planets have a mass fraction of water of greater 
than 40%. The terrestrial planets in our solar system are 
thought to have formed from water-devoid planetesimals. 
Hence the planetesimals must have formed either before the 
snow line moved inward or a fter the snow line mov e d out - 
ward past the Earth's orbit. I Oka. Nakamoto"fc Tdal (|201ll 1 
included ice opa city as well as the sil icate opacity and find 
similar results to lGaraud fc Lir] (|2007T ). They find that there 
is a deficit of solid mass at the later times and forma- 
tion of water-devoid planetesimals is impossible. It is im- 
probable for the planetesimals to form before the snow line 
moves inwards because the timescale involved is very short. 
iMachida fc Abel (|2010h investigated the possibility that the 
Earth formed from sublimating icy planetesimals after the 
snow line moved outwards. However, the extremely low ob- 
served water content on Earth appears to be an unlikely 
outcome, because of the competition between sublimation 
and collision of the planetesimals. In this paper we attempt 
to address the inconsistency between the theoretical models 
for the snow line evolution and the observations of the solar 
system. 

Protoplanetary discs may have a region of low turbu- 
lence at the disc mid-plane. The MRI drives turbu lence but 
may be suppressed by a low ionisation fraction l)Gammid 
ll996l ; lGammie fc Menoulll998f l. The inner parts of the disc 
are hot enough to be thermally ionised but farther away 
from the central star, cosmic rays are the dominant source 
of ionisation and these can only penetrate the surface layers. 
The mid-plane layer, the dead zone, has no turbulence. A 
variation of the turbulence with height from the mid-plane 
has been considered in some previous snow line models (e.g. 
iKretke fc Linll2007h . However, the mid-plane layer had a suf- 
ficiently large viscos ity that a steady disc was still achieved. 
IKretke fc Lml (|2010l ) consider the vertical structure of a lay- 
ered disc and find that the location of the snow line is a little 
more ambiguous with the vertical stratification. The upper 
layers of the disc have a cooler temperature and so the snow 
line is closer to the star in the surfaces than at the mid- 
plane. For an accretion rate of 10 _8 MQyr _1 they predict 
a snow li ne radius of ar o und 1 AU. This is similar to that 
found by iGaraud fc Linl l|2007f ) and lOka. Nakamoto"fc~ Jdi 
for the fully turbulent disc model. 

If there is no turbulence within the dead zone, the disc 
cannot be in a steady state and it is this possibility that we 
consider here. Material accumulates in the dead zone which 
can then become gravitationally unstable. While the infall 
accretion rate is high, this can lead to the gravo-magneto 
disc instability. The disc spends the majority of the time 
with a dead zone in a quiescent state with a low accretion 
rate on to the star. However, the extra turbulence driven 
by the gravitational instability in the dead zone increases 
the disc temperature until the MRI is triggered. This causes 
a large accretion outburst on t o the star that is thought 
to ex p lain FU Orionis outbursts dArmitage. Livio fc P ringle 



200 ll ; IZhu. Hartmann fc GammkT feoiO: Martin fc Lubowl 



20111 '). At later times, when the infall accretion rate drops, 



there is not sufficient material in the disc for the gravo- 
magneto in stability to operate but a dead zone may still 
be present |Martin et alj|2012ah . Given this rather different 



(from a steady-state disc) configuration and behaviour, in 
this paper, we consider the evolution of the snow line in a 
time dependent disc model that includes a dead zone. 



2 PROTOPLANETARY DISC MODEL 

The material in an accretion disc orbits the central mass, 
M, with a Keplerian velocity at radius R with angular ve- 
locity n = y^GM/R 3 (|Lvnden-Bell fc Pringldl 19741 ; iPringiel 
1981). We use a one dimensional layered disc model de- 
sc ribed in iMartin fc Lubowl (|201ll ) and further developed 
in lMartin et al. ( 2012al ) to evolve the total surface density, 
T,(R,t) and mid-plane temperature, T c (R,t). We take a so- 
lar mass star, M = 1 Mq , with a disc that extends from a 
radius of R = 5Rq up to R — 40 AU. Turbulence is driven 
by the MRI in some regions of the disc and we parametrise 
this viscosity with a lShakura fc Sunvaevl (1 19731) a parame- 
ter that we take to be 0.01 (e.g. iBrandenburg et alj|l995l : 
IStone et~al . 1996) a lthough there is still som e uncertainty in 
this value (e.g. King. Pringle fc Livid 12007). The disc con- 
sists of a turbulent surface layer of surface density, E m and 
temperature T m . Where it exists, the dead zone has surface 
density E g = E — E m . The extent of the dead zone in the 
disc is determined by a critical magnetic Reynolds number, 
ReM,crit(R, z). We use the analytical approximations for the 
surface d ensity in the turbulen t layer given in equations 26 
and 27 of IMartin et al.l |2012bl ). The innermost parts of the 
disc are externally he ated by a flux of radiation from the 
central star (see also IChiang fc Goldreichl fl997). We take 
the temperature of the star T star = 3000 K and its radius to 

be -Rstar = 2 Rq. 

We consider a mode l for the collapse of a molecu- 
lar cloud on to the disc l|Armitage. Livio fc Pringld I2OO1I ; 
IMartin et all l2012al ). Initially the accretion rate is 2 x 
10 -5 Mq yr -1 and this decays exponentially on a timescale 
of 10 5 yr. Material is added to the disc at a radius of 35 AU. 
The dead z one can becom e self gravitating when the Toomre 
parameter (|Too mrc 1964) is sufficiently small, Q < 2. We 
approximate the radius of the snow line to be where the 
te mperature of the d isc is T anow = 170 K. The disc models 
of lLecar et al.l (|2006h suggest that this is a good approxima- 
tion for varying disc mass and opacity. 

In Fig. [T] we show the evolution of the snow line in 
two model discs, one that is fully turbulent (or equivalently 
^?eM,crit = 0) and one that has a dead zone determined by 
a critical magnetic Reynolds number of iieM.crit = 5 x 10 4 . 
The fully turbulent model shows a monotonic decrease of the 
snow line radius in time so that at late times the snow line 
is within the radius of the Earth's orbit at 1 AU. We do not 
obtain the increase o f the snow line radius a t smal ler accre- 
tion rates shown by lOka. Nakamoto fc Idal l|201ll ) because 
there is sufficient mass in the disc that the model remains 
optically thick at the radius where the snow line occurs. For 
example, at a time of 2 x 10 6 yr, the infall accretion rate 
on to the disc has dropped to less than 10 -13 Mq yr~ x but 
the accretion on to the star is 6 x 10~ 9 Mq yr -1 . The disc 
does not remain in a quasi-steady state as the accretion rate 
decreases as assumed in previous models. 

For the model with 7?eM,crit = 5 x 10 4 yr, initially there 
is no dead zone and the disc evolution proceeds as for the 
case with J?eM,crit = 0. As the infall accretion rate drops, 
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Figure 1. The evolution of the icy regions of a protoplanetary disc around a solar mass star. Left: A fully turbulent disc (i?eM,crit = 0). 
Right: A disc with a dead zone determined by a critical magnetic Reynolds number of -ReM.crit = 5 X 10 4 . Time begins at t = 10 5 yr. 
The shaded areas show the icy regions of the disc. Where a dead zone exists there is an inner icy region (pale shaded areas) as well as 
the outer icy region (dark shaded areas). In the innermost parts of the disc, the radiation flux from the central star heats the disc above 
the snow line temperature, T snow . The dotted line shows the current radius of the Earth's orbit. 



the disc cools and a dead zone forms in the protoplanetary 
disc. The inner parts of the disc remain hot because they 
are heated by the radiation from the central star. Within 
the dead zone, an icy region exists. The outer parts of the 
dead zone become self gravitating and this drives a small 
amount of turbulence that heats the disc above the snow line 
temperature (see the sketch of the two disc models in Fig. [5]). 
The very brief periods of time during which there is no inner 
icy region are when the disc is in outburst. Self gravity heats 
the disc sufficiently to trigger the MRI within the dead zone. 
The temperature of the disc increases significantly and there 
is an FU Orionis type outburst where the accretion rate on 
to the star increase s by several orders of magnitude (see 
IMartin et ai1l2012al ). 

As with the fully turbulent disc, the outermost parts 
of the disc are cool and icy. At late times the outer icy 
region reaches a minimum radius of about 3.1 AU, roughly 
consistent with current observations. In our solar system, the 
Earth and other terrestrial planets likely formed after the 
inner icy region passed inside their radius because they are 
water devoid. This means that the Earth must have formed 
after a time of about 10 6 yr. At this time, there is little infall 
accretion on to the disc but there is a significant amount 
of mass in the disc, a total of about 0.3 Mg. The surface 
density at R — 1 AU is high at 7.5 x 10 4 gcm~ 2 . This is in 
agreement with the standard minimum mass solar ne bular 
at this radius of about 1.7 x 10 3 g em" 2 dHavashilll98lr i. It is 
comparable to the model of lDeschl (|2007r ) based on the 'Nice' 
model of planet formation that predicts a surface density of 
5x 10 4 gem -2 at the radius of the Earth's orbit. At this time 
there was still sufficient mass in the disc for planetesimal and 
Earth mass planet formation from the water-devoid disc. 



3 DISCUSSION 

Our simplified model clearly contains a number of uncer- 
tainties. Note that because we took the upper value for 
Tsnow, our results show the upper limit of the extent of 
the icy regions. Also, we have assumed, as did previous 
researchers, that the radius of the snow line in the gas 
disc is equivalent to t he radius of the snow line of the 
planetesimal disc (e.g. lOka. Nakamoto"fc ldall201ll ). How- 
ever, once the planetesimals grow large enough, they may 
decouple from the gas disc, and follow quasi-circular or- 
bits. If they were to move over the snow line, then they 
would sublimate or condense at the water-ice evaporation 
front (jStevenson fc Lunind[l988h . The p rocess of sublima- 
tion c an occur on a short ti mescale (e.g. ISack fc Baragiolal 
1 19931 : iBrown Sz Bolinal 2007) bu t cond ensation is less well 
understood. Marseille Sz Cazauxl (|201ll ) find there is an ex- 
tended region beyond the snow line where both icy and 
bare grains coexist, for a radial distance of about 0.4 au. 
This should be investigated further in future work with a 
disc model tha t track s the evolution of the water. While 
ICiesla fc Cuzzil (|2006l ) investigated this phenomenon in a 
fully turbulent disc, such calculations should be extended 
to a disc that contains a dead zone. 

There are other unknown parameters in the model 
such as t he value of th e critical magnetic Reynolds num- 
ber (e.g. iFleminel l2000h . and the viscosity parameter a. 
IMartin et al.l (|2012al )" showed that in order to reproduce 
FU Orionis outbursts, i?eM,crit must be a few 10 4 . There 
is also some uncertainty about whether the surface den- 
sity in the turbulent surface layer is determined by Ohmic 
resistivity (as we have assumed) or if ambipolar diffu- 
sion or the Hall effect could play an important role (e.g. 
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A: Steady, fully turbulent disc 




B: Non-steady disc with a dead zone 




Figure 2. Sketches of the disc structure for the two models. Fig- 
ure A represents a steady disc that is fully turbulent. The tem- 
perature decreases monotonically with radius and only the outer 
regions of the disc are icy. Figure B represents a disc with a dead 
zone. The outer parts of the dead zone are self-gravitating, more 
massive and thus have a higher temperature.. There are two icy 
regions in this model. The outer icy region is farther out than that 
predicted by the fully turbulent model. The inner icy region exists 
within the dead zone. The shaded regions correspond to the inner 
(pale shaded) and outer (dark shaded) icy regions as described in 

Fig. m 



IPerez-Becker fc Chiand l201ll: iBai fc Stond [ioTll l . On this 



last point, however, 



Wardlc & Salmeronl 



find the 



Ohmic resistivity term provides an average value for the 
active layer surface density for a range of vertical magnetic 
fields. We should also note that there appears to be some in- 
consistency between the accretion rates predicted by these 
dead zone models and observed T Tauri rates. The work 
presented here is not intended to be a comprehensive study 
of the evolution of the snow line. Rather, the purpose is to 
show that the formation a dead zone in a protoplanetary 
disc can keep the outer snow line farther away from the 
central star. Thus, a model that includes a dead zone can 
provide a solution to the problem that the Earth must form 
from water-devoid planetesimals that may not be present at 
1 AU in a fully turbulent disc model. This should be inves- 
tigated further in future work with more detailed numerical 
simulations once some of the current issues have been re- 
solved. 

On a more speculative note, the inner icy region, not 
predicted by the turbulent disc model, could allow for the 
formation of icy planets or gas giants in the inner regions of 
exo-solar systems. For example, at a time of 6 x 10 s yr there 
is a mass of 6.7 Mj in the inner icy region of the disc. The 
cores of giant planets are usually assumed to form in the cool 
region beyond the outer snow line. For example, hot Jupiters 
are giant planets that are very close to their central star. 



They are thought to have migrated through the d isc to their 
current position (e.g. iMasset fe Papaloizoull2003r ). However, 
with the proposed disc model (including a dead zone) it is 
possible that they could form much closer to the star, and 
closer to their observed locations. The newly formed core 
would not migrate inwar ds very fast because it wou ld be 
within a dead zone (e.g. iMatsumura fe p"ud ritz 2007j) and 
the dead zone extends closer to the star than the icy region 
(see Fig.[5J. It is possible that the inner icy region could exist 
for around 10 7 yr leaving perhaps sufficient time for giant 
planet core formation. However, we note that in order to 
form icy planetesimals we need not only a low temperature, 
but also water vapour to be present. For example, if the cold 
gas is dry, then icy planetesimals will neither form, nor will 
ice be added to pre-existing planetesimals. The speculative 
possibility of the formation of hot Jupiters at their current 
location needs further investigation. 



4 CONCLUSIONS 

We have computed the evolution of the icy regions of a pro- 
toplanetary disc that contains a dead zone. The disc evolves 
from times of high infall accretion through FU Orionis out- 
bursts to late times when the infall accretion rate is negli- 
gible and planetesimals are thought to form. In a fully tur- 
bulent disc, the unique snow line passes inside the radius 
of the Earth's orbit causing problems for the formation of 
our water-devoid planet Earth. However, when a dead zone 
forms in the disc, it prevents mass flow through the disc 
and the outer regions become self-gravitating. The turbu- 
lence driven by self-gravity increases the temperature of the 
outer parts of the dead zone and thus the outer icy region 
is much farther from the star as shown in the sketches in 
Fig. [2] Our model also predicts an inner icy region during 
early times of disc evolution, while a dead zone is present. 
This moves inside the radius of the Earth's orbit after a time 
of around 10 6 yr. The outer icy region exists at all times but 
has a minimum radius of around 3.1 AU, allowing for the 
formation of our water-devoid planet at its current radius. 

The key point of the present work is that the inclu- 
sion of a dead zone in time-dependent protoplanetary disc 
models can significantly change the evolution of the snow 
line. We have shown that it could resolve the apparent con- 
tradiction of the formation of the Earth from water-devoid 
planetesimals and it also introduces the possibility for icy 
planet formation close to the central star. 
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